Theory of chaotic atomic transport in an optical lattice 
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A semiclassical theory of chaotic atomic transport in a one-dimensional nondissipative optical 
lattice is developed. Using the basic equations of motion for the Bloch and translational atomic 
variables, we derive a stochastic map for the synchronized component of the atomic dipole moment 
that determines the center-of-mass motion. We find the analytical relations between the atomic 
and lattice parameters under which atoms typically alternate between flying through the lattice 
and being trapped in the wells of the optical potential. We use the stochastic map to derive 
formulas for the probability density functions (PDFs) for the flight and trapping events. Statistical 
properties of chaotic atomic transport strongly depend on the relations between the atomic and 
lattice parameters. We show that there is a good quantitative agreement between the analytical 
PDFs and those computed with the stochastic map and the basic equations of motion for different 
ranges of the parameters. Typical flight and trapping PDFs are shown to be broad distributions 
with power law "heads" with the slope —1.5 and exponential "tails". The lengths of the power 
law and exponential parts of the PDFs depend on the values of the parameters and can be varied 
continuously. We find analytical conditions, under which deterministic atomic transport has fractal 
properties, and explain a hierarchical structure of the dynamical fractals. 

PACS numbers: 42.50.Vk, 05.45.Mt, 05.45.Xt 



I. INTRODUCTION 

The transport properties of cold atoms in optical 
lattices depend on the lattice and atomic parameters 
and can be very diverse. The optical lattice is a peri- 
odic structure of micron-sized potential wells which is 
created by a laser standing wave made of two counter- 
propagating laser beams. Dilute atomic samples with 
negligibly small atom-atom interactions are used to 
probe experimentally single-atom phenomena. The 
atomic motion can take form of ballistic transport, 
oscillations in wells of the optical potential, Brown- 
ian motion, random walks, Levy flights, and chaotic 
transport. 

Most of the experimental 0, [3, i, i, 1, S and the- 
oretical P, @, H, m EH EH works on the atomic trans- 
port in optical lattices have been done in the context 
of laser cooling of atoms. One particularly interesting 
aspect of these works is the discovery of anomalous 
transport properties of cold atoms and Levy flights 
[T2, EH. A Levy flight is a random process in the space 
(time) domain with the distribution of the length (du- 
ration) of flight events that is given by a Levy law pos- 
sessing an algebraically decaying "tail" and infinite 
variance. As a consequence, the atomic trajectories 
(time series) have self-similar (fractal) nature and the 
probability to find superlong flights is not negligibly 
small as in a case of normal diffusion. 

In the context of sub-recoil laser cooling, the Levy 
flights have been found experimentally [lj, Q in the 
distributions of trapping and escape times for ultra- 
cold atoms trapped in a momentum state close to the 
dark state. It was shown in Refs. p], that not only 
the variance but the mean time for atoms to leave the 



trap is infinite. 

The anomalous atomic transport of cold atoms in 
optical lattices is of a different nature. A cold atom in 
an optical lattice can be trapped in the potential wells 
and it can move over many wavelengths in dependence 
on whether its energy is below or above the potential 
barrier. Transitions between these events are stochas- 
tically induced by the cooling force and its fluctuations 
(heating due to randomness of spontaneous emission 
and fluctuations of the atomic dipole moment). Mea- 
surement of the trajectory of a single cold ion in a 
one-dimensional optical lattice, by tracing its position 
through fluorescent photons, have demonstrated Levy 
flights Q. By deer casing the optical potential depth, 
a change of the transport characteristics from diffusive 
to quasiballistic was observed [1| . 

The anomalous properties of atomic transport dis- 
cussed briefly above are associated with random 
atomic recoils in spontaneous emission processes 
which are generally inevitable in any cooling scheme 
and which make the transport looks like a random 
walk. However, the problem may be considered not 
in a laser cooling but in a more general context as a 
study of deterministic motion of atoms with compar- 
atively large momentum interacting with a standing- 
wave light field. It is a fundamental nonlinear interac- 
tion between different degrees of freedom (the transla- 
tional and internal atomic degrees of freedom and the 
field ones) that can be treated in a Hamiltonian form 
when neglecting any losses. In this context sponta- 
neous emission may be considered as a noise imposed 
on the coherent atomic dynamics. It has been pre- 
dicted in Refs. [HI, EH that, besides the well-known 
transport properties of atoms in optical lattices, there 
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should exist a deterministic chaotic transport with a 
complicated alternation of atomic oscillations in wells 
of the optical potential and atomic flights over many 
potential wells when the atom may change the direc- 
tion of motion many times. This phenomenon looks 
like a random walk but it should be stressed that it 
may occur without any random fluctuations of the lat- 
tice parameters and any noise like spontaneous emis- 
sion. The deterministic chaotic transport is a result 
of chaotic atomic dynamics in the standing-wave field 
that means an exponential sensitivity of the inter- 
nal and translational atomic variables to small vari- 
ations in the lattice parameters and/or initial con- 
ditions. The chaotic transport and its manifestations 
like atomic dynamical fractals may occur both in clas- 
sical [II, [13, EGl] and quantized light fields [HEH- 
Spontaneous emission events interrupt the coherent 
atomic dynamics in random instants of time and may 
give rise anomalous statistical properties of atomic 
transport [22T |. 



In our previous papers [Tfl, UJ&, [2l|, l23| we have 
found chaotic atomic transport and dynamical frac- 
tals in numerical experiments and studied its proper- 
ties and manifestations. The ranges of the lattice and 
atomic parameters and initial conditions, for which 
the center-of-mass motion may be chaotic, have been 
established. In this paper we develop a semiclassical 
Hamiltonian theory of the chaotic atomic transport in 
a one-dimensional optical lattice and confirm the ana- 
lytical results by the numerical simulation. In Sec. II 
we derive the basic equations of motion and give the 
result of computation of the maximum Lyapunov ex- 
ponent whose positive values determine the ranges of 
the atom-field detuning A and initial atomic momen- 
tum po for which chaotic transport occurs. Sec. Ill 
briefly reviews distinct regimes of motion. Using ap- 
proximate solutions of the basic equations, we con- 
struct in Sec. IV a stochastic map for the synchro- 
nized component of the atomic dipole moment u that 
determines the chaotic transport. In Sec. V we intro- 
duce a simple illustrative model of random walking of 
the quantity arcsinit on a circle. Depending on the 
relations between the lattice and atomic parameters, 
the transport properties may be diverse. We find in 
Sec. V the condition under which the probability den- 
sity functions (PDFs) for the flights and trappings of 
the atoms either purely exponential or have promi- 
nent power law slopes. We derive the PDFs analyti- 
cally and compare the results with simulation of the 
stochastic map and the basic equations. The results 
obtained in Sec. IV are used in Sec. VI to find the con- 
ditions for appearing dynamical atomic fractals and to 
explain their structure. Finally, Sec. Vllgivcs conclu- 
sions. 



II. HAMILTON-SCHRODINGER 
EQUATIONS OF MOTION 

We consider a two-level atom with mass m a and 
transition frequency uj a , moving with the momentum 
P along the axis X in a one-dimensional classical 
standing laser wave with the frequency ujf and the 
wave vector kf. In the frame, rotating with the fre- 
quency ujf, the Hamiltonian is the following: 

P 2 1 

H = 1 — h(ui a — (jJf)cr z — Kl (<r_ + <7 + ) cos kfX. 

(i) 

Here <t±, z are the Pauli operators which describe the 
transitions between lower, |1), and upper, |2), atomic 
states, n is the Rabi frequency which is proportional 
to the square root of the number of photons in the 
wave y/n. The laser wave is assumed to be strong 
enough (n ^> 1), so we can treat the field classically. 
The simple wavefunction for the electronic degree of 
freedom is 



|*(t)) =o(t)|2)+6(t)|l), 



(2) 



where a and b are the complex-valued probability am- 
plitudes to find the atom in the states \2) and |1), 
respectively. Using the Hamiltonian fl}, we get the 
Schrodingcr equation 

da 
i — 
dt 

db 



Clb cos ktX, 



(3) 



—b — Via coskfX 



dt 2 

Let us introduce instead of the complex-valued prob- 
ability amplitudes a and b the following real-valued 
variables: 

u = 2Re (ab*) : v = -2 Im (ab*) , 

2 2 ( 4 ) 

z^\a\ 2 -\b\ 2 , 

where u and v are a synchronized (with the laser field) 
and a quadrature components of the atomic electric 
dipole moment, respectively, and z is the atomic pop- 
ulation inversion. 

In the process of emitting and absorbing photons, 
atoms not only change their internal electronic states 
but their external translational states change as well 
due to the photon recoil. If the atomic mean mo- 
mentum is large as compared to the photon momen- 
tum hkf, one can describe the translational degree of 
freedom classically. The position and momentum of a 
point-like atom satisfy classical Hamilton equations of 
motion. Full dynamics in the absence of any losses is 
now governed by the Hamilton- Schrodingcr equations 
for the real-valued atomic variables 

X = W r P, 

p = — it sin x, 

u = Av, (5) 
v = —Au + 2z cos x, 
z = —2v cos x. 
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where x = kfX and p = P/hkf are classical atomic 
center-of-mass position and momentum, respectively. 
Dot denotes differentiation with respect to the dimcn- 
sionless time t = fit. The normalized recoil frequency, 
u> r = hk 2 j/m a Q, <C 1, and the atom-field detuning, 
A = (ujf — oj a )/il, are the control parameters. The 
system has two integrals of motion, namely the total 
energy 



the plot marks the value of the maximum Lyapunov 
exponent A. In white regions the values of A are al- 



UJ r 2 A 

ti = — p — u cos X z, 



and the Bloch vector 



z 2 = l. 



(0) 



(7) 



The conservation of the Bloch vector length immedi- 
ately follows from Eqs. (@|. 

Equations of motion similar to the set ([5]) were ob- 
tained in our previous papers fbil . [Tol Il7j in order to 
describe the interaction between a two-level atom and 
the cavity radiation field in the strong-coupling limit. 
Taking into account a back reaction of the atom on 
the radiation field and within the semiclassical ap- 
proximation, we were able to get the corresponding 
version of the Hamilton-Schrbdinger equations for the 
atomic position x, momentum p, population inversion 
z, and two combined atom-field variables which were 
denoted by the same letters u and v as the atomic 
dipolc-moment components in Eqs. (T5T> Those equa- 
tions has been shown in Refs. [TJ, [l5|, QjJ to have 
a positive Lyapunov exponent in a wide range of the 
control parameters and initial atomic momentum pq. 
It implies dynamical chaos in the usual sense of ex- 
ponential sensitivity to small changes in initial condi- 
tions and/or control parameters. The same should be 
valid with the set of equations (O describing the dif- 
ferent physical situation — a two level atom in open 
space with a strong standing-wave field. 

Equations constitute a nonlinear Hamiltonian 
autonomous system with two and half degrees of free- 
dom which, owing to two integrals of motion, move on 
a three-dimensional hypcrsurface with a given energy 
value H. In general, motion in a three-dimensional 
phase space in characterized by a positive Lyapunov 
exponent A, a negative exponent equal in magnitude 
to the positive one, and zero exponent. The sum of 
all L yap unov exponents of a Hamiltonian system is 
zero [24j . The maximum Lyapunov exponent charac- 
terizes the mean rate of the exponential divergence of 
initially close trajectories, 



A = lim A(t), 



A(r)= lim -In ^4, (8) 



and serves as a quantitative measure of dynamical 
chaos in the system. Here, d(r) is a distance (in the 
Euclidean sense) at time r between two trajectories 
close to each other at initial time moment r = 0. The 
result of computation of the maximum Lyapunov ex- 
ponent in dependence on the detuning A and the ini- 
tial atomic momentum pq is shown in Fig. [TJ Color in 
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FIG. 1: Maximum Lyapunov exponent A vs atom-field 
detuning A (in units of the laser Rabi frequency Q) and 
initial atomic momentum po (in units of the photon mo- 
mentum hkf): Lo r = 10~ 5 , ito = zq = 0.7071, vo — 0. 

most zero, and the atomic motion is regular in the 
corresponding ranges of A and po- In shadowed re- 
gions positive values of A imply unstable motion. 

In all numerical simulations we shall use the normal- 
ized value of the recoil frequency equal to uj r — 10 -5 . 
The initial atomic position is taken to be xq = 0. 
The detuning A will be varied in a wide range, and 
the Bloch variables are restricted by the length of the 
Bloch vector ([7]). It should be noted that we use in 
this paper the normalization to the laser Rabi fre- 
quency f2, not to the vacuum (or single-photon) Rabi 
frequency as it has been done in our previous papers 
0j EH, EH, G~3 ■ So the ranges of the normalized con- 
trol parameters, taken in this paper, differ from those 
in the cited papers. Figure Q] demonstrates that the 
center-of-mass motion becomes unstable if the dimen- 
sionless momentum exceeds the value po f=a 300 that 
corresponds (with our normalization) to the atomic 
velocity v a ~ 1 m/s of a cesium atom in the field 
with the wavelength close to the transition wavelength 
A a ~ 852 nm. 



III. REGIMES OF MOTION 

A. Regular atomic motion at exact atom-field 
resonance 

The case of exact resonance, A = 0, was considered 
in detail in Ref. [l5|, [23| . Now we briefly repeat the 
simple results for the sake of self-consistency At zero 
detuning, the variable u becomes a constant, u = uq, 
and the fast (u, u, z) and slow (x, p) variables are sep- 
arated allowing one to integrate exactly the reduced 
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equations of motion. The total energy is equal to 
J r 2 



tt U r _ 

Ha = Y P 



ucosx, 



(9) 



where u = Uq. The center-of-mass atomic motion in 
this spatially periodic potential of the standing wave 
is described by the simple nonlinear equation for a free 
physical pendulum 



x + uj r UQ sin x = 0, 



(10) 



and does not depend on evolution of the internal de- 
grees of freedom. 

The translational motion is trivial when uq is zero. 
The atom moves in one direction with a constant ve- 
locity, and the Rabi oscillations arc modulated by the 
standing wave. Equations ([9]) and (|10[) describe the 
atomic motion in the simple cosine potential uq cos x 
with three types of trajectories which arc possible in 
dependence on the value of the energy H: oscillator- 
like motion in a potential well if Hp < uq (atoms are 
trapped by the standing- wave field [251] ). motion along 
the scparatrix if Hq = uq, and ballistic-like motion if 
Hq > uq. Exact solutions of Eq. (1101) are easily found 
in terms of elliptic functions (see |l5l. [23j). 

As to internal atomic evolution, it depends on the 
translational degree of freedom since the strength of 
the atom-field coupling depends on the position of 
atom in a periodic standing wave. At A = 0, it is 
easy to find the exact solutions of Eqs. © 



maximum Lyapunov exponent A depends both on the 
parameters uj r and A, and on initial conditions of the 
system It is naturally to expect that off the res- 
onance atoms with comparatively small values of the 
initial momentum pa will be at once trapped in the 
first well of the optical potential, whereas those with 
large values of po will fly through. The question is 
what will happen with atoms, if their initial kinetic 
energy will be close to the maximum of the optical 
potential. Numerical experiments demonstrate that 
such atoms will wander in the optical lattice with al- 
ternating trappings in the wells of the optical potential 
and flights over its hills. The direction of the center- 
of-mass motion of wandering atoms may change in a 
chaotic way (in the sense of exponential sensitivity 
to small variations in initial conditions). A typical 
chaotic atomic trajectory is shown in Fig. [5] 




v(t) = ±y/l — u 2 cos I 2 J cosxdr' 



Xo 



z(r) = =F\/l — u 2 sin I 2 J cosxdr' + xo 



(11) 



where u = Uo, and cos [x(t)] is a given function of the 
translational variables only which can be found with 
the help of the exact solution for x [Til . [23[ . The sign 
of v is equal to that for the initial value zq and 



Xo = T arcsm 



(12) 



is an integration constant. The internal energy of the 
atom, z, and its quadrature dipole- moment compo- 
nent v could be considered as frequency-modulated 
signals with the instant frequency 2cos[x(r)] and the 
modulation frequency u) r p(r), but it is correct only 
if the maximum value of the first frequency is much 
greater than the value of the second one, i. e., for 
\u) r p \ < 2. 



B. Chaotic atomic transport off the resonance 

In Fig. [T] we show the A-map in the space of the 
initial momentum po and detuning A values. The 



FIG. 2: Typical atomic trajectory in the regime of chaotic 
transport: xo = 0, po = 300, zo = —1, uo — vq = 0, 
uj r = 10~°, A = —0.05. Atomic position is shown in units 
of kj 1 , time r in units of SI -1 . 

It follows from ([5]) that the translational motion 
of the atom at A ^ is described by the equation 
of a nonlinear physical pendulum with the frequency 
modulation 



x + LU r u(r) sinx = 0, 



(13) 



where u is a function of all the other dynamical vari- 
ables. 



IV. STOCHASTIC MAP FOR CHAOTIC 
ATOMIC TRANSPORT 

Chaotic atomic transport occurs even if the normal- 
ized detuning is very small, |A| -c 1 (Fig. [T]). Under 
this condition, we will derive in this section approx- 
imate equations for the center-of-mass motion. The 
atomic energy at |A| -C 1 is given with a good accu- 
racy by the simple resonant expression Returning 
to the basic set of the equations of motion wc may 
neglect the first term in the fourth equation since it 
is very small as compared with the second one there. 
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However, we cannot now exclude the third equation 
from the consideration. Using the solution (fTTjl for v, 
we can transform this equation as 



= ±Ay/l - u 2 



cos%, 



where 



2 / cos xdr' 



Xo- 



(14) 



(15) 



Far from the nodes of the standing wave, Eq. (fT4")l 
can be approximately integrated under the additional 
condition, \uj r p\ <C 1, which is valid for the ranges 
of the parameters and the initial atomic momentum 
where chaotic transport occurs. Assuming cos.t to 
be a slowly-varying function in comparison with the 
function cosx, we obtain far from the nodes the ap- 
proximate solution for the it-component of the atomic 
dipolc moment 



sin ± 



A 



2cosx 



sin x + C 



(16) 



where C is an integration constant. Therefore, the 
amplitude of oscillations of the quantity u for com- 
paratively slow atoms (|w r p| <C 1) is small and of the 
order of |A| far from the nodes. 

It follows from the third and forth equations in the 
set (O that u satisfies to the equation of motion for a 
driven harmonic oscillator with the natural frequency 
|A| and the driving force 2zcosx whose frequency is 
space and time dependent. At |A| = 0, the synchro- 
nized component of the atomic dipole moment u = is 
a constant whereas the other Bloch variables z and 
v oscillate in accordance with the solution (fTTj) . At 
| A. | 7^=0 and far from the nodes, the variable u per- 
forms shallow oscillations for the natural frequency 
|A| is small as compared with the Rabi frequency. 
However, the behavior of u is expected to be very 
special when an atom approaches to any node of the 
standing wave since near the node the oscillations of 
the atomic population inversion z slow down and the 
corresponding driving frequency becomes close to the 
resonance with the natural frequency. As a result, 
sudden "jumps" of the variable u arc expected to oc- 
cur near the nodes. This conjecture is supported by 
the numerical simulation. In Fig. [3] we show a typical 
behavior of the variable u for a comparatively slow and 
slightly detuned atom. The plot clearly demonstrates 
sudden "jumps" of u near the nodes of the standing 
wave and small oscillations between the nodes. 

Approximating the variable u between the nodes by 
constant values, we can construct a discrete mapping 
u m = f(u m —i), where u m is a value of u just after the 
TO-th node crossing (including multiple crossings of the 
same node). In order to estimate the values of u just 
after the crossing, one needs to integrate Eq. (fT4"|) near 
a node. Since u changes largely in a small region near 



U 0.6 




4000 



FIG. 3: Typical evolution of the atomic dipole-moment 
component u for a comparatively slow and slightly detuned 
atom: xo = 0, po = 550, v — 0, uo = zo = 0.7071, 
UJ r = 10~ 5 , A = -0.01. 



nodes (so small that the atomic momentum has no 
time to change its value significantly) we may use the 
Raman-Nath approximation of the constant velocity. 
In this approximation, we have 



X 



smx + xo- 



(17) 



Substituting this expression into Eq. (|14|) . we can 
integrate it in the interval x ^ tt (which comprises 
only one node) and obtain the value of u just after 
crossing the first node 



Ui ~ sin arcsinuo ± 



sin ^arcsin uq ± 
v J 



Att 



cos xdx 



(18) 



^rPnodc 



-VPnode 



zqEq 



^rPnodc 



where 



Pnode 



(19) 



is the value of the atomic momentum at the instant 
when the atom crosses a node (which is the same with 
a given value of the energy H for all the nodes), Jo and 
Eq are zero-order Bessel and Weber functions, respec- 
tively. In the limit of the large argument 2/ (u> r pnode), 
both the functions have a harmonic asymptotics [27j, 
and the expression (fl8|) reduces to the form 



u\ w sin ± 



A 



X Vq cos 



~ 4 



zq sin 



^rPnodo 4 

Zq] + arcsinuo) . 



(20) 
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It is the deterministic solution obtained with the ap- 
proximations mentioned above. It should be stressed 
that the solution contains trigonometric functions 
with large values of the arguments which arc in- 
versely proportional to the atomic momentum. In the 
Raman- Nath approximation, we take p ~ Pnodo- In 
fact, even small deviations from this mean value may 
result in large changes in the magnitude of the trigono- 
metric functions. Therefore, they can be treated as, 
practically random variables in the range [—1,1]. Be- 
yond the Raman-Nath approximation, the value of the 
atomic momentum p depends on the value of u m which 
changes every time when the atom crosses a node. So, 
we can replace arguments of the trigonometric func- 
tions by random variables. Finally, we introduce the 
stochastic map 




(21) 

where 4>m are random phases to be chosen in the range 
[0, 2ir]. When deriving this map, we neglected the 
term zq in Eq. (|20p which is small as compared with 
the factor yj ir/aj r p no< i e . 

With given values of A, u) r , and p no de, the map 
([2"Tj) has been shown numerically to give a satisfactory 
probabilistic distribution of magnitudes of changes in 
the variable u just after crossing the nodes. The 
stochastic map ([21]) is valid under the assumptions 
of small detunings (|A| <§; 1) and comparatively slow 
atoms (|w r p| <C 1). Furthermore, it is valid only for 
those ranges of the control parameters and initial con- 
ditions where the motion of the basic system ([5]) is un- 
stable. For example, in those ranges where all the Lya- 
punov exponents are zero, u becomes a quasi-periodic 
function and cannot be approximated by the map. 

The stochastic map (|2"Tj) allows to reduce the basic 
set of equations of motion ([5]) to the following effective 
equations of motion: 

x = uj r p, 

p = — u m smx, (22) 
to = Lu r p nodc |<5(cosx)|, 

where u m is found from Eq. (|2ip . The third equation 
in the set (|22|) gives a correspondence between the 
continuous evolution of the atomic motion and a dis- 
crete crossing number to. The integration of the delta 
function 5(cosx) over time at points with cosx = 
gives ±l/(w r pnode) in dependence on the direction of 
motion and whether the serial number of the node is 
even or odd. Since we calculate absolute values of the 
delta function, it is easy to show that to is a constant 
if cos a; ^ and increases by one if cosx = 0. 



V. STATISTICAL PROPERTIES OF 
CHAOTIC TRANSPORT 

A. Model for chaotic atomic transport 

With given values of the control parameters and the 
energy H, the center-of-mass motion is determined by 
the values of u m (see Eq. (fT3|)). One can obtain from 
the expression for the energy @ the conditions under 
which atoms continue to move in the same direction 
after crossing a node or change the direction of motion 
not reaching the nearest antinode. Moreover, as in the 
resonance case, there exist atomic trajectories along 
which atoms move to antinodes with the velocity going 
asymptotically to zero. It is a kind of scparatrix-likc 
motion with an infinite time of reaching the stationary 
points. 

The conditions for different regimes of motion de- 
pend on whether the crossing number to is even 
or odd. Motion in the same direction occurs at 
(— l) m+1 u m < H, scparatrix-likc motion - - at 
{-l) m+1 u m = H, and turns — at {-l) m+1 u m > H. 
It is so because even values of to correspond to cos x > 
0, whereas odd values — to cos a; < 0. The quantity u 
during the motion changes its values in a random-like 
manner (see Fig. [3]) taking the values which provide 
the atom either to prolong the motion in the same di- 
rection or to turn. Therefore, atoms may move chaot- 
ically in the optical lattice. The chaotic transport 
occurs if the atomic energy is in the range < H < 1. 
At H < 0, atoms cannot reach even the nearest node 
and oscillate in the first potential well in a regular 
manner (see Fig.[Tj). At H > 1, the values of u are al- 
ways satisfy to the flight condition. Since the atomic 
energy is positive in the regime of chaotic transport, 
the corresponding conditions can be summarized as 
follows: at |u| < H, atom always moves in the same 
direction, whereas at \u\ > H, atom either moves in 
the same direction, or turns depending on the sign of 
cos a; in a given interval of motion. In particular, if 
the modulus of u is larger for a long time then the 
energy value, then the atom oscillates in a potential 
well crossing two times each of two neighbor nodes in 
the cycle. 

The conditions stated above allow to find a di- 
rect correspondence between chaotic atomic transport 
in the optical lattice and stochastic dynamics of the 
Bloch variable u. It follows from Eq. (|2"Tj) that 
the jump magnitude u m — u m -i just after crossing 
the TO-th node depends nonlinearly on the previous 
value u m -i- For analyzing statistical properties of 
the chaotic atomic transport, it is more convenient to 
introduce the map for arcsinu m 

9 m = arcsinii m = 



= A, / smc/) m + arcsinu m _!, 

V ^rPnodc 

where the jump magnitude does not depend on a cur- 
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rent value of the variable. The map (f23|) visually looks 
as a random motion of the point along a circle of unit 
radius (Fig. 2]). The vertical projection of this point is 
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FIG. 4: Graphic representation for the maps of u m and 
9 m = arcsin u m . H is a given value of the atomic energy. 
Atoms either oscillate in optical potential wells (trapping) 
or fly through the optical lattice (flight). 

u m . The value of the energy H specifies four regions, 
two of which correspond to atomic oscillations in a 
well, and two other ones to ballistic motion in the 
optical lattice. 

We will call "a flight" such an event when atom 
passes, at least, two successive antinodes (and three 
nodes). The continuous flight length L > 2n is a dis- 
tance between two successive turning points at which 
the atom changes the sign of its velocity, and the dis- 
crete flight length is a number of nodes I the atom 
crossed. They are related in a simple way, L ~ nl, for 
sufficiently long flight. 

Center-of-mass oscillations in a well of the optical 
potential will be called "a trapping". At extremely 
small values of the detuning (the exact criterion will 
be given below), the jump magnitudes are small and 
the trapping occurs, largely, in the 2-7r-wide wells, i. e., 
in the space interval of the length 2tt. At intermediate 
values of the detuning, it occurs, largely, in the 7r-widc 
wells, i. c. in the space interval of the length tt. Far 
from the resonance, |A| > 1, trapping occurs only in 
the 7r-wide wells. Just like to the case of flights, the 
number of nodes I, atom crossed being trapped in a 
well, is a discrete measure of trapping. 



then the internal atomic variable 9 m = arcsin u m 
just after crossing the m-th node may take with the 
same probability practically any value from the range 
[— 7r/2,7r/2] (see Fig. 3}. With given values of the 
recoil frequency uj r = 10~ 5 and the energy in the 
range < H < 1 corresponding to chaotically moving 
atoms, large jumps take place at medium detunings 
|A| ~ 0.1. The probability P_ for an atom to turn 
just after crossing a node is equal to the probability 
to get to one of the trapping regions in Fig. 2] (which 
one depends on whether the crossing number m is even 
or odd) and is given by 

arccos -ff 1 
P- = < 5 , P+ = l-P-, (25) 

where P + is the probability to prolong the motion in 
the same direction after crossing the node. It is easily 
to get from Eq. ([25]) the probability for an atom to 
cross I successive nodes before turning 
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It is a flight probability density function (PDF) in 
terms of the discrete flight lengths. The exponential 
decay means that the atomic transport is normal for 
sufficiently large values of the jump magnitudes of the 
variable u. 

The statistics of the center-of-mass oscillations in 
the potential wells can be obtained analogously. With 
large values of the jump magnitudes (|24p . trapping 
occurs, largely, in the 7r-wide wells. The probability 
for a trapped atom to cross the corresponding well 
node I times before escaping from the well is 



P tr (l)=P_ l P + 



arccos H 
1 | exp 
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(27) 



C. Statistics of chaotic atomic transport at 
small jump magnitudes of u 

In this subsection we consider the case of small val- 
ues of the jump magnitudes of the variable u 
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B. Statistics of chaotic atomic transport at large 
jump magnitudes of u 

If the jump magnitudes of the variable u are suffi- 
ciently large 
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Now, it may take a long time for an atom to exit from 
one of the trapping or flight regions in Fig. [4] So, 
we need to calculate the time of exit of the random 
variable 6 m = arcsin u m from one of these regions. 
The result will depend on what is the length of the 
corresponding circular arc in Fig. [5] as compared with 
the jump lengths. 

Firstly, let us consider the case if the jump lengths 
are small as compared with the lengths both of the 
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flight and trapping arcs, i. e., if 



<C min{arcsin-ff, arccosff}. (29) 
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Motion of 9 m along the circle can be now treated as a 
one-dimensional diffusion process for a fictitious par- 
ticle described by the equation 



dV{9,m) _ D d 2 V(6,m) 



dm 



ae 2 



(30) 



where V is a probability density to find the particle 
at the angular position 9 just after crossing the m-th 
node and D is the corresponding diffusion coefficient 
of the particle 
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(31) 

which was calculated for the particle jumping ran- 
domly with the mean magnitude (9 m — m _i) (equal 
to zero in our case of symmetric jump distribution) 
and the variance ((9 m — # m -i) 2 )- Thus, we reduced 
the task to the first passage time probability problem 
for a continuous Markov process. Using the results 
of this theory for a Wiener diffusion process [28[ de- 
scribed by Eq. (f3"Tj) . we calculate the probability den- 
sity for a particle to exit from the interval 6 C ± # max 
after crossing I nodes 
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. (32) 

where 9q is an initial angular position of the particle in 
the region under consideration and 9 C is the center of 
the region (with four possible values 0, it/2, it, 3n/2). 

Since in a case of small jumps a particle gets to the 
region near its limit point, it is possible to replace 9q 
by the quantity (9o = 9 C — 9 max + e with a small positive 
value of e. Expanding the cosine in (|32p in the vicinity 
of the nodes and taking into account only the terms 
of the first order of smallncss, we get the flight and 
trapping PDFs 
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where Q is a normalization constant that is the same 
in both the cases. When deriving the first and second 



expressions (|33|) . we replaced max by arcsin H and 
arccos H, respectively. It is easily to realize that when 
j exceeds the value # max / (jry/Dl) , the corresponding 
terms of the series (|33| decrease rapidly. 

The PDFs (f3"3"| can be written in a much more sim- 
ple form for two conditions imposed on the number 
of jumps (node crossings) I a particle needs to quit 
the interval 9 C ± 9 max . If I > 9 2 aax /D, then all the 
terms in the sums (|33D arc small as compared with 
the first one. Both the flight and trapping statistics 
are exponential in this case. 

To the contrary, if I <C 9 2 aaK /D, then one should 
take into account a large number of terms in the sums 
(1331) , and each sum can be replaced approximately by 
the integral. In this case the result does not depend 
on the length of the region 2# max and we get the power 
law decay 
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both for the flight and trapping PDFs. The power- 
law statistics (|34[) implies anomalous atomic trans- 
port. Dividing the length 2# max by the jump magni- 
tude (|28p. we obtain the quantity 
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(35) 



that is a minimum number of jumps (node crossings) 
a particle needs to pass the region through. With the 
number of jumps I < l cl , a particle, randomly moving 
on the circle, may get out of the interval 9 C ± # max 
only through the same border where it got in. The 
statistics of exit times is known to be a power-law 
one with the transport exponent equal to —1.5. The 
length of the interval does not matter in this case. 
If I ^ 'cr ; then particles may exit through both the 
borders, and the corresponding statistics cannot be 
approximated by a single transport exponent. If / > 
l 2 r = 9 2 aax /D 1 then we expect an exponential decay. 

The size of the trapping and flight regions depends 
on the value of the atomic energy H (see Fig. @|. At 
H > V2/2 (arcsin H > tt/4), the flight PDF has a 
longer decay than the trapping PDF. On the contrary, 
at H < \/2/2, the ftr's decay is longer than the iVs 
one. If the jump magnitude (|28[) is of the order of the 
size of the flight or trapping regions 
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arcsin H <C — , 



arccos H <C — , 



(36) 



then a particle may pass through the region making 
a small number of jumps I. So, the approximation of 
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the diffusion process (|29|) fails, and the corresponding 
PDF is exponential. 
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trapping events. They are normalized histograms for 
a number of standing-wave nodes I which the atom 
crossed in a run. Each PDF is computed with a single 
but very long (up to t ~ 10 8 for the basic equations 
of motion) atomic trajectory. 
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FIG. 5: The flight P fl and trapping P tr PDFs for a chaot- 
ically moving atom with comparatively small values of 
the jump magnitudes of the variable u. The detuning 
A = —0.001 is small, and the energy value H — 0.724 
(po = 535) provides approximately equal sizes of the flight 
and trapping regions in Fig. [4] White and black circles rep- 
resent results of integration of the basic ((5]l and reduced 
(|22l) equations of motion, respectively, and the solid lines 
represent the analytical PDFs (|33[) . %o = 0, zq = — 1, 
Wo = = 0, to,- — 10~ 5 . 

In order to check the analytical results obtained in 
this section, we compare them with numerical sim- 
ulation of the reduced (|22|) and basic (O equations 
of motion. The reduced equations (f2"2"| describe the 
centcr-of-mass motion modulated by a jump-like be- 
havior of the internal variable u which satisfies to the 
stochastic map (f2"Tj) . To explore different regimes of 
the chaotic atomic transport we integrate Eqs. (|22"f 
and ([5]) numerically with two values of the detuning 
A = —0.001 and A = —0.01 (the cases of small and 
medium jump magnitudes of the variable u, respec- 
tively) and different values of the atomic energy and 
momentum and compute the PDFs of the flight and 



FIG. 6: The same as in Fig. [5] but with the energy 
H = 0.101 (po = 402) providing a narrow flight region 
as compared with the trapping one. 

In Fig. [5]we compare the results (in a log-log scale) 
in the case of small jump magnitudes of the variable 
u (A = —0.001) and approximately equal sizes of 
the flight and trapping regions (H = 0.724 ~ \/2/2). 
White and black circles represent the results of numer- 
ical integration of the basic (O and reduced (|2"2"j) equa- 
tions of motion, respectively, and the solid line repre- 
sents the analytical predictions ([55)) . The agreement 
between the PDFs computed with Eqs. ©, (|2"2"|) . and 
is rather good. All the flight and trapping PDFs 
exhibit in their dependence on the crossing number 
/ three kinds of decay which are fairly approximated 
by the formulas (|33[) . The critical value of the cross- 
ing number I is estimated to be l cr ~ 55 with chosen 
values of the parameters and initial conditions. In 
the range / < 55 the PDFs are expected to demon- 
strate the power law decay with the slope —1.5 given 
by the formula flM|)- In the range 55 < I < l 2 CI ~ 3000, 
there should exist a number of different transport ex- 
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ponents. The initial part of this range is fitted by the 
same power-law function. However, the other part of 
the range cannot be fitted by a simple function. In 
the range I > 1% T — 3000, the decay is expected to 
be purely exponential in accordance with the the first 
expression in Eqs. (f3"3")) . 

In Fig. [6] we show the PDFs computed at the same 
detuning A = —0.001 as in Fig. [5] but with the value 
of the atomic energy H = 0.101 for which the flight 
region in Fig. U is narrower than the trapping one. 
The jump magnitude is now of the order of the length 
of the flight region (see (J2U)) and the flight PDFs 
are expected to be exponential (l cr ~ 9). The crit- 
ical value of the crossing number I for trapping is 
Z cr ~ 58 with given values of the energy and the de- 
tuning. Therefore, in the range I < 58 the trapping 
PDFs are expected to vary as the power law decay 
with the slope —1.5 given by the formula ([3"4")) . In 
the range I > 1% T — 3400, the decay is expected to 
be purely exponential in accordance with the the sec- 
ond expression in Eqs. (|53"|) . The trapping PDFs in 
the range 58 < I < 1% X ~ 3400 are neither power law 
nor exponential demonstrating a complicated behav- 
ior. This speculation is confirmed by the numerical 
computation shown in Fig. [6l 

In order to demonstrate what happens with larger 
values of the jump magnitudes, we take the detun- 
ing to be A — —0.01 increasing the jump magnitude 
in ten times as compared with the preceding cases. 
With the taken value of the energy H = 0.8055 we 
provide a slight domination of flights over trappings. 
The jump magnitude is now so large that particles 
may pass both through the flight and trapping re- 
gions making a small number of jumps. It is expected 
that all the PDFs, both the flight and trapping ones, 
should be practically exponential in the whole range 
of the crossing number I. It is really the case (see 
Fig.©. 



VI. DYNAMICAL FRACTALS 

Various fractal-like structures may arise in chaotic 
Hamiltonian systems [III [29|, [3(j . In our previous pa- 
pers 

nam 

we have found numerically fractal proper- 
ties of chaotic atomic transport in cavities and optical 
lattices. In this section we apply the analytical re- 
sults of the theory of chaotic transport, developed in 
the preceding sections, to find the conditions under 
which the dynamical fractals may arise. 

We place atoms one by one at the point xq = 
with a fixed positive value of the momentum po and 
compute the time T when they cross one of the nodes 
at x = — 7r/2 or x = 3ir/2. In these numerical experi- 
ments we change the value of the atom-field detuning 
A only. All the initial conditions po — 200, zq = —1, 
= Vo = and the recoil frequency uj r = 10 -5 are 
fixed. The exit time function T(A) in Fig. [S] demon- 
strates an intermittency of smooth curves and compli- 
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FIG. 7: The flight P fl and trapping P tr PDFs for a chaot- 
ically moving atom with comparatively large values of 
the jump magnitudes of the variable u. The detuning 
A = —0.01 is medium, and the energy value H = 0.8055 
(po = 550) provides a domination of the flight events over 
the trapping ones. White and black circles represent re- 
sults of integration of the basic ((5} and reduced (|22[) equa- 
tions of motion, respectively, and the solid lines represent 
the analytical PDFs (|33[) . £o = 0, zo = — 1, uq = vq = 0, 

LU r = 10~ 5 . 



cated structures that cannot be resolved in principle, 
no matter how large the magnification factor. The 
second and third panels in Fig. [5] demonstrate succes- 
sive magnifications of the detuning intervals shown 
in the upper panel. Further magnifications reveal 
a self-similar fractal-like structure that is typical for 
Hamiltonian systems with chaotic scattering [26|, [3(| ■ 
The exit time T, corresponding to both smooth and 
unresolved A intervals, increases with increasing the 
magnification factor. Theoretically, there exist atoms 
never crossing the border nodes at x = —tt/2 or 
x = 3tt/2 in spite of the fact that they have no obvious 
energy restrictions to do that. Tiny interplay between 
chaotic external and internal atomic dynamics pre- 
vents those atoms from leaving the small space. The 
similar phenomenon in chaotic scattering is known as 
dynamical trapping. In Ref. [l7| we have computed 
the Hausdorff dimension for the similar fractal and 
shown that it is not an integer. We note that the 
fractal-like structures similar to that shown in Fig. [5] 
arise in numerical experiments with longer distances 
between the border nodes but the corresponding com- 
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FIG. 8: Fractal-like dependence of the time of exit of 
atoms T from a small region in the optical lattice on the 
detuning A: po — 200, zo = — 1, = «o = 0. Magnifica- 
tions of the detuning intervals are shown. 



putation time increases significantly with increasing 
this distance. 

Various kinds of atomic trajectories can be char- 
acterized by the number of times to atom crosses 
the central node at i = ir/2 between the border 
nodes. There are also special separatrix-like trajec- 
tories along which atoms asymptotically reach the 
points with the maximum of the potential energy, hav- 
ing no more kinetic energy to overcome it. In differ- 
ence from the separatrix motion in the resonant sys- 
tem (A = 0) with the initial atomic momentum p cr , 
a detuned atom can asymptotically reach one of the 
stationary points even if it was trapped for a while in 
a well. We specify the mS-trajectory as a trajectory 
of an atom crossing the central node m times before 
starting the separatrix-like motion. Such an asymp- 
totical motion takes an infinite time, so the atom will 
never reach the border nodes. 

The smooth A intervals in the first-order structure 
(Fig. [8j upper panel) correspond to atoms which never 
change the direction of motion (m = 1) and reach the 
border node at x = 3tt/2. The singular points in the 
first-order structure with T = oo, which are located 
at the border between the smooth and unresolved A 
intervals, are generated by the lS'-trajectories. Anal- 
ogously, the smooth and unresolved A intervals in the 
second-order structure (second panel in Fig. [SJ corre- 



spond to the 2-nd order (m = 2) and the other trajec- 
tories, respectively, with singular points between them 
corresponding to the 2«S'-trajectories, and so on. 

There are two different mechanisms of generation of 
infinite exit times, namely, dynamical trapping with 
an infinite number of oscillations (to = oo) inside the 
interval [— tt/2, 3tt/2] and the separatrix-like motion 
(m ^ oo). The set of all the values of the detunings, 
generating the separatrix-like trajectories, was shown 
to be a countable fractal in Rcfs. [ij], HH, whereas 
the set of the values generating dynamically trapped 
atoms with m = oo seems to be uncountable. The 
exit time T depends in a complicated way not only 
on the values of the control par ameters but on initial 
conditions as well. In Ref. |23| we presented a two- 
dimensional image of this function in two coordinates 
of the initial atomic momentum pq and the atom-field 
detuning A. 

The fractal-like structure with smooth and unre- 
solved components may appear if atoms have an al- 
ternative either to turn back or to prolong the motion 
in the same direction just after crossing the node at 
x = tt/2. For the first-order structure in the upper 
panel in Fig. [51 it means that the internal variable u of 
an atom, just after crossing the node for the first time 
(cos a; < 0), satisfies either to the condition U\ < H 
(atom moves in the same direction), or to the con- 
dition u\ > H (atom turns back). If u\ = H, then 
the exit time T is infinite. The jumps of the vari- 
able u after crossing the node are deterministic but 
sensitively dependent on the values of the control pa- 
rameters and initial conditions. We have used this fact 
when introducing the stochastic map in Sec. IV. Small 
variations in these values lead to oscillations of the 
quantity arcsinui around the initial value arcsinuo 
with the amplitude \A\y/ir/ (w r Pnodo ) ■ If this ampli- 
tude is large enough, then the sign of the quantity 
Ui — H alternates and we obtain alternating smooth 
(atoms reach the border x = 3tt/2 without chang- 
ing their direction of motion) and unresolved (atoms 
turns a number of times before exit) components of 
the fractal-like structure. The singular points at the 
border of such components correspond to the atoms 
with infinite exit time. 

If the values of the parameters admit large jump 
magnitudes of the variable u (see condition (|24p). then 
the dynamical fractal arises in the energy range < 
H < 1, i. e., at the same condition under which atoms 
move in the optical lattice in a chaotic way. In a 
case of small jump magnitudes, fractals may arise if 
the initial value of an atom uq is close enough to the 
value of the energy H, i.e., the atom has a possibility 
to overcome the value u = H in one jump. Therefore, 
the condition for appearing in the fractal T(A) the 
first-order structure with singularities is the following: 



| arcsinuo ~~ arcsinifj < |A 



(37) 
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The formation of the second-order structure is ex- 



12 



plained analogously. If an atom made a turn after 
crossing the node for the first time, then it will cross 
the node for the second time. After that, the atom 
either will turn or cross the border node at x = —w/2. 
What will happen depend on the value of U2- How- 
ever, in difference from the case with m = 1, the con- 
dition for appearing an infinite exit time with m = 2 
is U2 = — H. Furthermore, the previous value u\ is 
not fixed (in difference from u ) but depends on the 
value of the detuning A. In any case we have u\ > H 
since the second-order structure consists of the trajec- 
tories of those atoms which turned after the first node 
crossing. Near the singular borders of the second- 
order structure we have u\ ~ H, whereas far from 
them u\ may be much larger. In order for an atom 
would be able to turn after the second node crossing, 
the magnitude of its variable u should change suffi- 
ciently to be in the range < —H. The atoms, 
whose variables u could not "jump" so far, leave the 
space [— 7r/2, 3tt/2]. The singularities are absent in the 
middle part of the second-order structure shown in the 
second panel in Fig. [5] because all the corresponding 
atoms left the space after the second node crossing. 
The variable U2 oscillates with varying A generating 
oscillations of the exit time. In vicinity of the borders 
of the second-order structure the required jump mag- 
nitude U2 — U\ is smaller and is equal approximately 
to 2H . Some atoms could turn for the second time or 
follow separatrix-like trajectories which generate sin- 
gularities of the exit-time function. So, the condition 
for appearing singularities in the second-order struc- 
ture is the following: 



2arcsini? < |AL/ — - — . (38) 
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With the values of the parameters taken in the simu- 
lation, we get the energy H = 0.2 + A/2. It is easy to 
obtain from the inequality (|38[) the approximate value 
of the detuning |A| ss 0.0107 for which the second- 
order singularities may appear. In the lower panel in 
Fig.[S]one can see this effect. No additional conditions 
are required for appearing the structures of the third 
and the next orders. 

It should be noted that the inequality (|38p is op- 
posite to the first part of the inequality (|2"5)) that de- 
termines the condition for appearing power law de- 
cays in the flight PDF. Therefore, dynamical fractal 
may appear in those ranges of the control parameters 
where the Levy flights are impossible and vice versa. 
However, the trapping PDF may have a power law 
decay. It should be noted that the inequality (|38f in 
difference from (|37p is strongly related with the chosen 
concrete scheme for computing exit times. It is not re- 
quired with other schemes, say, with three antinodes 
between the border nodes. 



VII. CONCLUSION AND DISCUSSION 

We have studied scmiclassical transport of atoms 
in a one-dimensional optical lattice in the framework 
of the idealized model of the atom-field interaction. 
The center-of-mass motion is strongly affected by the 
evolution of the atomic Bloch variables — the syn- 
chronized and quadrature components of the atomic 
dipole moment u and v and the population inversion 
z. There are ranges of the atom- field detuning A and 
the initial momentum po where the atomic transport 
has been shown to be chaotic with a complicated al- 
ternation of trappings in wells of the optical potential 
and flights over them. We developed a semiclassical 
theory of the chaotic atomic transport in optical lat- 
tices in terms of a random walk of the Bloch variable 
u and confirmed the analytical results by direct nu- 
merical simulation of the basic and reduced equations 
of motion. 

Based on a jump-like behavior of the internal vari- 
able u for atoms crossing nodes of the standing laser 
wave, we constructed a stochastic map that deter- 
mines the center-of-mass motion. We found the re- 
lations between the detuning A, recoil frequency u> r , 
and the atomic energy H, under which atoms may 
move in the lattice in a chaotic way. To illustrate 
statistical properties of the chaotic transport, we pro- 
posed a simple model of random walking of the quan- 
tity arcsinw on a circle. 

The statistical properties of chaotic atomic trans- 
port strongly depend on the relation between the jump 
magnitude of the variable u and the atomic energy. 
Both the length of flight and the duration of trapping 
are measured by a number of nodes I atom crossed 
during the flight and trapping events, respectively. 
The PDFs for the flight P A (l) and trapping P tI (l) 
events were analytically derived to be exponential in a 
case of large jumps. In a case of small jumps, the kind 
of the statistics depends on additional conditions im- 
posed on the atomic and lattice parameters, and the 
distributions Pr(1) and Pt r (0 were analytically shown 
to be either practically exponential or functions with 
long power law parts with the slope —1.5 but exponen- 
tial "tails". The comparison of the PDFs computed 
with analytical formulas, the stochastic map, and the 
basic equations of motion has shown a good agreement 
in different ranges of the atomic and lattice parame- 
ters. We used the results obtained to find the ana- 
lytical conditions, under which the fractal properties 
of the chaotic atomic transport were observed, and to 
explain the structure of the corresponding dynamical 
fractals. 

Since the period and amplitude of the optical po- 
tential and the atom-field detuning can be modified 
in a controlled way, the transport exponents of the 
flight and trapping distributions are not fixed but 
can be varied continuously, allowing to explore differ- 
ent regimes of the atomic transport. Our analytical 
and numerical results with the idealized system have 
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shown that deterministic atomic transport in an op- 
tical lattice cannot be just classified as normal and 
anomalous one. We have found that the flight and 
trapping PDFs may have long algebraically decaying 
parts and a short exponential "tail" . It means that in 
some ranges of the atomic and lattice parameters nu- 
merical experiments reveal anomalous transport with 
Levy flights (see Figs. [5] and [6]). The transport expo- 
nent equal to —1.5 means that the first, second, and 
the other statistical moments arc infinite for a reason- 
ably long time. The corresponding atomic trajecto- 
ries computed for this time are self-similar and frac- 
tal. The total distance, that the atom travels for the 
time when the flight PDF decays algebraically, is dom- 
inated by a single flight. However, the asymptotical 
behavior is close to normal transport. In other ranges 
of the atomic and lattice parameters, the transport is 
practically normal both for short and long times (see 
Fig. [3). 

In conclusion we should emphasize that the theory 
of chaotic atomic transport in an optical lattice pre- 
sented in this paper does not take into account spon- 
taneous emission which is inevitable in any real exper- 
iment. Speculation about small changes in the atomic 
momentum of the order of hkf (caused by sponta- 
neous emission) as compared with its initial values of 
the order of 300 -f- 500hkf used in our computations 



is not a convincing argument because the atomic mo- 
mentum is equal to zero at turning points, and sponta- 
neous emission may change the center-of-mass motion 
in vicinity of these points. Spontaneous emission can 
be incorporated in the equations of motion of the type 
of our basic equations (O [13] by the quantum Monte 
Carlo method [3l[. In general, we expect a competi- 
tion between the coherent but chaotic atomic dynam- 
ics and random processes of spontaneous emission. In 
particular, the stochastic map (|2"Tj) constructed in this 
paper should be generalized to include jumps of the 
variable u to the zero value at random time moments. 
What will happen with the flight and trapping PDFs 
is an open question. We plan to study the atomic 
transport in an optical lattice with spontaneous emis- 
sion in the future. 
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